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ABSTRACT 

Kinematical parameterisations of disc galaxies, employing emission line observations, 
are indispensable tools for studying the formation and evolution of galaxies. Future 
large-scale Hi surveys will resolve the discs of many thousands of galaxies, allowing 
a statistical analysis of their disc and halo kinematics, mass distribution and dark 
matter content. 

Here we present an automated procedure which hts tilted-ring models to Hi data 
cubes of individual, well-resolved galaxies. The method builds on the 3D Tilted Ring 
Fitting Cbde (TiRiFiC) and is called FAT (Fully Automated TiRiFiC). 

To assess the accuracy of the code we apply it to a set of 52 artificial galaxies 
and 25 real galaxies from the Local Volume Hi Survey (LVHIS). Using LVHIS data, 
we compare our 3D modelling to the 2D modelling methods DiskFit and ROTCUR. 

A conservative result is that FAT accurately models the kinematics and the 
morphologies of galaxies with an extent of eight beams across the major axis in the 
inclination range 20°-90° without the need for priors such as disc inclination. When 
comparing to 2D methods we find that velocity fields cannot be used to determine 
inclinations in galaxies that are marginally resolved. We conclude that with the 
current code tilted-ring models can be produced in a fully automated fashion. This 
will be essential for future Hi surveys, with the Square Kilometre Array and its 
pathhnders, which will allow us to model the gas kinematics of many thousands of 
well-resolved galaxies. Performance studies of FAT close to our conservative limits, 
as well as the introduction of more parameterised models will open up the possibility 
to study even less resolved galaxies. 

Key words: galaxies: ISM, galaxies: kinematics and dynamics, galaxies: structure, 
methods: data analysis, surveys 


1 INTRODUCTION 

The accurate description of the kinematics of galaxies is cru¬ 
cial to get insight into how they form and evolve during their 
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lifetimes. The Doppler-shifted Hi line is one of the main trac¬ 
ers of these kinematics; Hi is ubiquitous, largely unaffected 
by absorption and often extends far further out than other 
probes. As such, parameterisation of the Hi distribution and 
dynamics has been done for several decades, mostly on an 
individual galaxy basis and in a highly interactive fashion. 
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Successful and objective kinematic parametrisation can 
lead to a wealth of understanding in galaxy formation and 
evolution. For exam ple, it can be used to investigate the 
origin of warps (e.g. Jozsa l2007l : Ivan der Kruit fc FreemanI 


I 2 OIII : iHaan &: Braun 20141 ') and with that, gas accretion 


via in trinsic and large-scale spin alignment lIPoDPing et al.l 
l2014h . but also the study of dark an d baryonic matter dis¬ 
tributions (ISwaters et al.l[2012l . l20l3) . 

Current Hi surveys such as the We sterbork Hi Survey o f 
Irregular and Spiral galaxies (WHISP) (Ivan der Hulst 2002), 


the Local Volume Hi Survey (LVHIS) (Koribal^ 
DiskMass (lMartinssonll201ll), ATLAS^° (Serra et al 


(lMartinssonll20L^ 

lue disc Project (|\ 


201011 

I2OI2II 


and the Blue disc Project Iwang et al. I I2OI3II have already 
observed tens to hundreds of galaxies for which kinematic 
parametrisation can be obtained. Although galaxy models 
can s till be produced int eractively for this number of sys¬ 
tems JSwaters et al.ll2012h . the need for automation is obvi¬ 
ous. 

This need is further emphasised when looking at the 
planned H i surveys for the Au stralian SKA Pathfinder 
(ASKAP) lljohnsto n et al.| I2OO8II and the APERture Tile 
In Focus (Apertif) ( Oosterloo et al.l [ 2009 I I . The Wide field 
ASKAP L-band Legacy All-sky Blind surveY (WALLABY; 
iKoribal^ I2OI2II and the Westerbork Northern Sky Hi 
Survey (WNSHS) will provide an all-sky Hi survey at ~ 
15” — 30” resolution, a significant improvement upon pre¬ 
vious single-dish surveys. WALLABY and WNSHS are ex¬ 
pected to resolve thousands of galaxies to an extent that 
their kinematical parameters can be retrieved llDuffv et al.l 
l2012li . This paper is the first in a series that presents an al¬ 
gorithm for automatically extracting kinematic parameters 
for these detections. 

Owing to the fact that in general the orbits of the neu¬ 
tral gas in disc galaxies deviate only slightly from circu¬ 
larity, the common approach to parameterise this geome- 
try and kinematics i s to make use of the tilted-ring model 
dRogstad et al.lll974h . In this framework, gas discs are mod¬ 
elled as a sequence of rings of increasing radius, where each 
ring has its own geometry and kinematics. Historically, the 
fitting was performed on velocity fields constructed from 
the Hi emission line observations. When analysing velocity 
fields, the tilted-ring model is typically under-constrained, 
such that the fitting process needs to be guided interactively 
(introducing, often subjective, priors to effectively regularise 
the fitting process). In a manifold of Hi studies a process is 
followed where the velo city field of a galaxy is fitted with 
the G IPSY task rotcur llBegemarJI 19871 : van der Hulst et al.l 
Il992^ in several iterations. A typical fit starts by varying sev¬ 
eral parameters simultaneously to obtain initial guesses, and 
these parameters are subsequ ently fitted one at a time until 
they converge to a solution (iBegemanI Il987l : Ide Blok et al.l 

l2008h . 

This approach has several problems. As the velocity 
field is used for constructing the model, beam smearing and 
projection effects can severely affect the final outcome. Al¬ 
ready early on it was recognised that the inner rotation curve 
slope is biased low, due to beam smearing, when the data 
have le ss than 7 res olution elements inside the Holmberg 
radius (iBosmal Il978f) . At high inclinations the projection 
effects become so severe that reliable rot ation curves can 
no longer be retrieved fro m velocity fields dSancisi fc AllenI 
ll979l:ISofue fc Rubirj200lli . Even though this effect has been 


long known, an accurate investigation into the exact inclina¬ 
tion at which projection effects become limiting to accurate 
analysis has not been undertaken. 

The coupling between the rotational velocities and in¬ 
clination can also affect tilted ring fits of velocity fields. As 
we only observe velocities along the line-of-sight, the ob¬ 
served velocities relate to the rotational velocities as Vobs = 
Wrot X sin(i). In particular, inclination and rotation veloc¬ 
ity become mathematically degenerate in velocity fields of 
galaxies with solid-body rotation curves, as well as at low 
inclinations where the line-of-sight component of the rota¬ 
tion is similar in amplitude to the disc velocity dispersion. 
As a result, galaxies with inclinations i < 40°are typically 
not well-s uited for tilted r ing mo delling employing velo city 
fields fe.g. iBegemanlUoSTl . but see iBershadv et aLil201Cll ~) . 

In recent years multiple codes have become available 
that are optimised for different aspects of 2D kinematic anal¬ 
ysis. For example, the tilted-ring model has been adapted for 
integral field unit (IFU) observatio ns of elliptical galaxie s 
at optical wavelengths (kinemetry. iKrainovic et al. l2()()6l ~) . 
This technique deals with the underconstraint by first ap¬ 
plying a grid based approach to find the global minima 
and subsequently applies a minimisation. Another tech¬ 
nique is to break the degeneracy outright by fitting flat 
disc models to the data, an approach adopted by DiskFitQ 
llSnekkens fc Sellwoodl[2007h . 

Hi data as obtained with modern interferometers are 
mostly 3D, with two spatial axes and one spectroscopic 
axis. However, when constructing velocity fields the obser¬ 
vations are compressed into a 2D representation of the ve¬ 
locities which can cause a loss of information. Addition¬ 
ally, the method of co mpression, can add another bias to 
the data analysis (e.g. Ide Blok et al.l[2008l') . By fitting the 
3D observations directly several problems, such as beam 
smearing and projection effects, can be overcome. It is also 
possible to fit instances where the disc crosses the line of 
sight multiple times. The most extreme cases of the lat¬ 
ter example are edge-on galaxies where the disc crosses 
the line-of-sight an infinite number of times. While rota- 
tion curves alone may be extracted usin g other methods 
llSofuelllsi^ : iKregel fc van der Kruitll2004l '). 3D fitting tech¬ 
niques are most suitable for extracting kinematical informa¬ 
tion from such galaxies. 

3D studies can not only negate projection effects but 
they can also be used to better under stand the structure 
of galaxies fe.g. lVoigtlander et al.ll2oi3 l. The disadvantages 
of 3D fitting are that the fitting process is computationally 
expensive an d more sensitive t o inhomogeneities in the gas 
distribution d Jozsa et al.ll200^ L A model describing the full 
3D data requires additional parameters which in turn need 
to be determined or fixed, therefore the models may still be 
under-constrained and an interactive fitting procedure has 
so far been required. _ 

The Tilted Ring Fitting Code (TiRiFiC IJozsa et al.l 
1200711^ has long been the only published and freely available 
3D tilted-ring fitting code. The current version of TiRlFiC 
has an unprecedented degree of flexibility, owing to the fact 
that past development has been focused on enabling detailed 


^ http://www.physics.rutgers.edu/-spekkens/diskfit/ 
^ http://gigjozsa.github.io/tirific/ 
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multi-layer modelling. Models employing a large range of pa¬ 
rameters are required for the interpretation of ultra-deep Hi 
observations of nearby galaxies, such as those in the Hy- 
drogen Accre tion in LOcal GAlaxieS (HALO GAS) Survey 
jHeald^£_^ 201l|_ Zschaechner et ^ I2OIII: iGentile et al.l 


l2013l : lKamphuis et al.ll2013l : lde Blok et al.ll201'j ). Due to the 

complexity of these gas distributions and the wide range of 
TiRlFiC model parameters allowed, TiRlFiC requires inter¬ 
active fitting and little emphasis has been put on developing 
a truly automated fitting method. Here, we present our re¬ 
cent development to overcome this weakness of TiRlFiG in 
the limit where all model parameters are directly determined 
from the data and no priors are used. 

In this paper, we present and test an IDL code that 
utilises the strengths of TiRlFiC without the need for in¬ 
teraction during the fitting processes. We test the code on 
a set of 52 artificial galaxies and 25 real Hi observations se¬ 
lected from LVHIS. In addition we explore the differences 
between the 3D and 2D approaches and different methods 
of stabilising the fits by comparing the 3D fits with rotcur 
and DiskFit models. FAT does not utilise all of TiRiFiC’s 
flexibility but aims to reliably fit an asymmetric warp, sur¬ 
face brightness profile and a symmetric rotation curve to a 
range of different galaxies. As such its output can be used to 
retrieve these kinematical parameters for large surveys or it 
can be used as an objective starting point for more complex 
models. This code is a prototype for the 3D branch of the 
WALLABY kinematical pipeline. 

This paper is structured as follows. In § [5] we detail 
the philosophy for this pipeline and specify the initial input 
and structure of the code presented in this paper. § [3] intro¬ 
duces the test modules we have created for the code and § 
13] presents the results of these tests. Finally we will discuss 
these results and suggest possible improvements to the code 
in §[S]and conclude and summarise in §[3| 


2 A KINEMATIC MODELLING PIPELINE 

Figure [T] shows a conceptual flow diagram for the pipeline 
that is being developed to model the kinematics of resolved 
detections from WALLABY and WNSHS. The key principle 
of this pipeline is that it should automatically extract kine¬ 
matical parameters in a uniform and reliable manner. This 
paper will focus on the 3D branch of the pipeline, whereas 
the 2D branch will be presented in an accompanying paper 
llOh et al jliH). 

The input for the pipeline can be either a velocity field 
or a data cube, depending on the initial estimates that are 
extracted from the data through the source finder. Based on 
these estimates the pipeline will perform (more efficient) 2D 
fitting in cases where this approach is thought to be reliable 
and (more direct) 3D fitting for more complicated cases. If 
the 2D fit fails the pipeline will try to fit the full data cube 
in the 3D branch. 

The 3D branch presented in this paper is Fully Auto¬ 
mated TiRlFiC ( FAT) which is cu rrently an IDL wrapper 
around TiRiFiC (Ijozsa et al.ll2007h . The aim of FAT is to 
extract from an Hi line observation the surface brightness 
profile (SBR), position angle (PA), inclination (INCL) and 
rotational velocity (Vrot) as a function of radius and the 
scale height (ZO), intrinsic dispersion (SDIS) and central co¬ 



Figure 1. Conceptual flow diagram of the complete WALLABY 
kinematic modelling pipeline. 


ordinates (RA, DEC, Vgya) as single values without any hu¬ 
man interaction. In the following we describe the pipeline 
inputs and the fitting procedure. 

Due to the continuously improving computational capa¬ 
bilities FAT is designed with the philosophy that accurate 
results are far more important than speed. Even so, in this 
section we first investigate the speed bottlenecks in a sys¬ 
tematic way. In TiRlFiC there are two major ways in which 
the fitting process can be significantly slowed. The first is 
having a large amount of voxels in the data cube. The other 
bottleneck is related to how TiRlFiC creates the models. 
A model is created by populating the model geometry with 
Hi clouds, i.e. point sources, in a Monte-Carlo fashion and 
subsequen tly convolving thi s point source model with a 3D- 
Gaussian lljozsa et al.ll2007li . The user can control the num¬ 
ber of point sources used in the model with a specific input 
parameter that assigns an approximate flux to each point 
source. A large number of point sources, i.e a low flux for 
each point source, will result in a large increase in computing 
time. However, a low number of point sources will increase 
the discretisation noise and results in lower accuracy. 

In order to determine the optimum speed settings with¬ 
out loss of accuracy we have tested a single TiRlFiC itera¬ 
tion on five artificial galaxies. In this test TiRlFiC performs 
a minimisation through sequential fi tting of all param - 
eters using the golden-section method ijPress et al.l ll99l L 
The artificial galaxies were created with TiRlFiC and white 
noise with random amplitude was added to mimic observ¬ 
ing uncertainties. In these fits we varied the number of pixels 
across the minor axis FWHM of the beam from 1-4 as well 
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Figure 2. Comparison of speed vs. accuracy. From left to right 
panels show duration, i and their average as determined from a 
simple fit of five galaxies created with TiRlFiC on a grid of pixels 
per minor axis beam and number of point sources per model. The 
images are smoothed for display purposes. 


as the number of point sources used to generate the model. 
TiRlFiC fits each parameter sequentially until all parame¬ 
ters are fitted in a so-called loop. It then evaluates whether 
specific convergence criteria are satisfied. If not, it will start 
again from the first parameter. We timed the duration of 
every fit and the average duration of a single TiRlFiC loop 
and averaged these two quantities. The reason for consider¬ 
ing both is that in TiRlFiC the number of loops required for 
convergence can change with changing accuracy. This rea¬ 
soning was confirmed by small differences in the time grids 
build up from a single loop and a full iteration. However, as 
this depends on the complexity of the fit/galaxy we consider 
both. The fit and loop times are normalised for each galaxy 
before averaging. 

Figure [2] shows the grid for, from left to right, nor¬ 
malised duration, normalised final of the model and a 
speed-accuracy grid in which the former two are averaged. 
The grid consists of a 4 (pixels per beam) by 10 (number of 
point sources) roster. Formally the optimal settings for this 
set of galaxies, i.e. the minimum of the combined duration 
X^ grid, is at 4 pixels per beam and 1x10® point sources 
per model. However, it is clear that accurate results with 
3 pixels across the minor axis FWHM of the beam require 
at least 3x 10® point sources per model. Additionally we can 
see from the left most panel in Figure[^that for models with 
more than 2x 10® point sources, the fitting slows down signif¬ 
icantly without a large increase in accuracy. FAT therefore 
adjusts the point source flux such that the number of point 
sources per model is always between 5x10® and 2.2x10®. 
Regridding the input cubes is optional but in this paper we 
always use data with 4 pixels across the minor axis FWHM. 


2.1 FAT Input 

FAT requires an input catalogue identifying a set of galaxies 
to be fitted, a configuration file and a data cube following 


the FITS standard for each galaxy. All fitting is done in ob¬ 
served units, e.g. arcsec not kpc, and thus a distance is not 
required. However, if a galaxy distance is provided outputs 
are converted to physical units such as kpc and Mq. Pro¬ 
viding a distance can also help with the initial estimates, 
especially the initial scale heights and scale lengths as they 
are based on physical empirical estimates. 

TiRlFiC requires initial estimates for the kinematic pa¬ 
rameters to be fitted. The initial estimates for the veloc¬ 
ity width, the size and galaxy centre are obtained from the 
source finder SoFiA0 JSerra et al.ll2015h . SoFiA also pro¬ 
vides a mask that is used for the creation of an integrated 
moment map and intensity weighted mean velocity field, i.e. 
moment 0 and moment 1 maps. The reason for choosing 
SoFlA is twofold; it provides uniform initial estimates and 
masks from which to start the fitting process and it is a 
product of the WALLABY source finding group, thus en¬ 
suring compatibility of the pipeline with the WALLABY 
source finder. 

The moment maps created from the SoFlA mask are 
used to estimate the PA. This estimate is the average of 
the PA obtained from a simple ellipse fit to the integrated 
moment map and the PA determined from the angle be¬ 
tween the maximum and minimum velocity in the intensity 
weighted mean velocity field. Following the estimate of the 
PA, the inclination is determined as the average of the axis 
ratio at three different fractions of the peak intensity (i.e. 
at i X Ipeak, I X Ipeak, ^ X Ipeak) and five different PAs as 
determined from the error on the PA estimate, i.e PA-error, 
PA-error/2., PA, PA-|-error/2 and PA-|-error. As mentioned 
above, these preparatory steps are fully integrated into the 
FAT pipeline. 


2.2 The 3D Fitting Procedure 

FAT is a wrapper around SoFlA and TiRlFiC that enables 
fully automated tilted-ring fits. TiRlFiC is run a number of 
times with increasing model complexity. In addition, FAT 
adds the following features to the TiRlFiC functionality: i) 
estimating the disk’s outer radius, ii) regularising (smooth¬ 
ing) the parametrisation, iii) assessing fit quality/model con¬ 
vergence, iv) controlling fitting speed, v) estimating errors. 
In the following we describe these steps in detail. While 
TiRlFiC is being consta ntly developed, m ost features used 
in FAT are described in ijozsa et alJ (l2007l L In addition we 
utilise TiRlFiC’s ability to fit two independent half disks 
to an observation. The TiRlFiC stopping criteria are deter¬ 
mined by the SATDELT parameter, which varies through¬ 
out the different steps in the code. In the code presented 
here TiRlFiC is always used in golden section mode. 

The first step in the fitting process (see Figure jS]), af¬ 
ter initial estimates are obtained, is to fit a fiat disc, i.e. a 
model in which all parameters except the brightness and the 
rotational velocity of each ring have the same value at all 
radii. We adopt the major axis FWHM of the beam as the 
ring width in order to produce a stable fit and to ensure that 
neighbouring rings are independent. The first ring however 
is merely one fifth of a FWHM in order to keep the central 


® https://github.com/SoFiA-Admin/SoFlA/ 
This paper uses vO.4.0 
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Figure 3. Flow chart for FAT. 


rotational velocity at zero without introducing artefacts. Ad¬ 
ditionally we halve the ring size for very small galaxies, i.e. 
less than six beams across the major axis, and double the 
ring size for large galaxies, i.e. more than 50 beams across 
the major axis, beyond the inner 20 resolution elements. 

In this first step the main aim is to get accurate central 
coordinates for the model. For the other parameters we aim 
to ht values that place them within the local minima of their 
optimal values. We look for convergence on the central po¬ 
sition, such that between two subsequent iterations it does 
not change by more than the size of a voxel. If convergence 
is not achieved the code will rerun the fitting process with 
the updated central coordinates. Additionally, FAT checks 
that the other parameters remain within the boundaries of 
the initial estimates and updates them when the initial es¬ 
timates appear incorrect (see Figure [3]). If the code fails to 
converge within 50 iterations, the fitting procedure is halted 
and flagged as failed. 

After a satisfactory ht for a hat disc is obtained, FAT 
proceeds to ht radial variations for the inclination, PA, ro¬ 
tational velocities and the surface brightness. In addition to 
introducing the radial variations we split the model along 
the minor axis and ht both sides independently except for 
the rotation curve. For small galaxies (< 6 beams in diam¬ 
eter) we do not vary the inclination and PA radially. 

In this step the code evaluates several criteria before 


accepting the ht. Firstly, the internal TiRlFiC convergence 
criteria must be satished instead of TiRlFiC hnishing be¬ 
cause the maximum of ten loops through all parameters was 
reached. If the latter is the case the code will average the ra¬ 
dially varying values and try again. Secondly, there must be 
no more than two consecutive rings close to the boundaries 
determined in the initial estimates. If there are the bound¬ 
aries are expanded and TiRlFiC is run again. Finally the 
SBR of the last ring must be in an acceptable range. This is 
determined by comparing the SBR in the last ring to a c utoff 
value derived using Equation 4 from Ijozsa et all ll2007ll : 


CTTiRiFiC 


Jy km s 


= 9.0 X 10" 


Hr 


Ilbe 


( 1 ) 


where orimFiC is the cutoff value for each ring, Ilbeam 
is the beam solid angle and firing is the solid angle of 
the rin£|. This formula is independent of the noise value 
and therefore we scale it by the ratio of the noise in the 
fitted cube to tha t in the test cubes of 3.6 mjy beam~^ 
djozsa et al.l 1200^ 1. We consider the model to be converged 
when the surface brightness of the last ring is between 1-5 
times this cutoff value. If it is less than the cutoff value the 
last ring is dropped, whereas if it is more than 5 times the 
cutoff value the model is extended by one ring. Whenever 
these iterations lead to a number of rings that has been 
fitted before, for example by adding a ring and subtracting 
a ring in the following iteration, or differs by more than 2 
rings from the initial estimate we fix the number of rings to 
remain at this number. 

Once the previously stated criteria are satished we 
smooth the parameters with a polynomial of order 0, 2, 3, 4 
or 5. This is because the simple ring to ring htting routine 
of TiRlFiC can cause the parameters to vary artihcially 
from ring to ring in a saw tooth pattern and this leaves 
the possibility for large outliers. The 1st order is excluded 
because a linear increase from centre to outer radius is not 
a physical representation of PA and INCL. To determine 
the order of the polynomial to use for each parameter we 
ht, through a Monte Carlo approach, all polynomials and 
adopt the one with the lowest reduced value. After the 
smoothing the parameters are rehtted, with the ring to ring 
variations kept constant, and once more the brightness of 
the last ring is compared to the cutoff value. 

As a hrst empirical attempt at estimating errors on 
tilted ring models, FAT calculates errors for the htted 
parameters on a ring by ring basis. These errors are 
estimated by the maximum of either the variation in the 
Monte Carlo hts or the difference between the smoothed 
and unsmoothed values. We want to stress that this is a first 
empirical attempt at determining errors and is by no means 
a formal or statistically correct error. Part of the investi¬ 
gation here is to see whether such errors perform reasonably. 


^ Note that in Ijozsa et al.l ||2007^ the exponent appears to be 
solely applied to the ring solid angle however this is an error. 
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3 TEST MODULES 

3.1 Artificial Galaxies 

We have constructed 52 artificial galaxies to test the accu¬ 
racy of FAT. We build the models around two different ro¬ 
tation curves (see Figure |3|: one resembling a slowly rising 
rotation curve typical of intermediate-mass spiral galaxies 
and one base d on a smoothed version of the rotation curve 
of NGC 891 (lOosterloo et al.|l2007l : iFraternali et aLll201ll l 
in order to represent massive spiral galaxies. For each type 
of galaxy we adopt an exponential surface brightness pro¬ 
file with a scale length of 10 kpc, which then gets locally 
perturbed to mimic surface brightness variations. This pro¬ 
file is truncated at a specific radius by decreasing the scale 
length of the profile by a factor of 20. The values of other 
parameters of the base galaxies, such as inclination and size, 
were set to optimal values for tilted ring fitting. The final 
model observations are then derived from these base models 
by varying one parameter. As we have set the base galaxy 
values optimally, the effect of the parameter that is var¬ 
ied can be isolated. The values for the base galaxies are 
shown in columns two and three of Table [T] and the values 
to which they were varied are in the fourth column. We cre¬ 
ate these artificial galaxies as a distribution of point sources, 
add white noise and subsequently smooth them to produce 
a range in the number of beams across the major axis. 

In these models the FWHM of the beam is kept con¬ 
stant at 30” and the physical truncation radius is kept at 25 
and 20 kpc for the intermediate and massive galaxy, respec¬ 
tively. In order to vary the number of beams that span the 
diameter defined by the truncation radius the distance is ad¬ 
justed. Afterwards this distance is used to convert the scale 
heights from fixed physical values to angular scale heights. 

The noise is calculated as a fraction of the average sig¬ 
nal in the noise-free cube; first in each channel the average 
is calculated of all pixels that exceed 10% of the maximum 
signal in the cube. The average of all these averages is sub¬ 
sequently taken as the average signal in the cube. 

In order to warp some of the test galaxies we vary the 
angular momentum vector of the initial disc linearly with 
radius, in a random direction. Table [T] lists, in the third col¬ 
umn, the variations of 9, the angle between the line-of-sight 
and the angular momentum vector, and (j), the angle between 
the North-South axis of the sky and the angular momentum 
vector. We keep the inner disc flat and start applying this 
variation to the angular momentum vector at | of the disc 
extent such that the maximum variation, as stated in Ta¬ 
ble [T] is reached at the cutoff radius. We then project the 
warped rings back to an inclination and PA for each ring. 
This leads to realistic warps with changes in PA and incli¬ 
nation of ~10°. 

In total, we produce 25 models for the high mass and 
intermediate mass galaxy. To study the effect of a reduced 
velocity resolution (in terms of resolution elements) when 
fitting dwarf galaxies, we add to the 50 models two further 
models resembling a dwarf galaxy (see Figure ^ with the 
parameters of the base galaxies. The exception to this is the 
surface brightness profile, which is produced anew with a 
scale length of 7.5 kpc and is truncated at a physical radius 
of 15 kpc. 


3.2 LVHIS Test Module 


Real galaxies can contain many deviations from highly sym¬ 
metrized artificial galaxies. We therefore run FAT on real 
observational data, comparing it to the 2D fitting methods 
ROTCUR and DiskFit. Since we do not know the correct so¬ 
lution a priori, the comparison mainly aims at finding sys¬ 
tematic differences between 2D and 3D methods. ROTCUR 
and TiRlFiC fit the same model in different dimensions and 
both codes need to be guided by hand to an acceptable re¬ 
sult or, as is attempted here, embedded in a wrapper that 
makes assumptions about the properties of real galaxies and 
applies them automatically. 

For DiskFit the problem is not under-constrained as 
it assumes the disc to be flat; this allows the code to fit 
global minima and to calculate realistic error bars. DiskFit 
can include streaming motions into the final model to sim¬ 
ulate bars or radial motion and also includes a model for a 
weak, quadratic, symmetric warp. Therefore, the compari¬ 
son between rotcur and DiskFit can be used to investigate 
systematic errors introduced by complex warps in the outer 
parts and bars in the inner parts of these galaxies. 

In order to test the pipeline’s 2D and 3D components we 
selected 25 galaxies from th e southern sample of the LVHIS 
catalogue dKoribalskil l201Cli) . LVHIS is a complete sample 
of galaxies selecte d from HIPASS llKoribalski et al] l2004l : 
iMever et al.ll20M') with a systemic velocity ^ 550 km s“^ 
(-^ 10 Mpc). The southern sample was further restricted to 
have 5 < —30°. This sample was observed with the Australia 
Telescope Compact Array (ATCA) with similar angular and 
velocity resolution to the proposed values for WALLABY 
(FWHM -30", channel width -4 km s’^). 

The 25 galaxies were selected by taking the LVHIS Hi 
mass function and selecting the first two most regularly ro¬ 
tating galaxies, as determined from a visual inspection of the 
velocity fields, from 0.1 dex-wide mass bins. This sample was 
interactively fitted w ith ROTCUR (for a standard description 
of the approach see lOh et al.l [200^1 . In particular, it was 
decided by inspection whether or not the galaxy kinemat¬ 
ics were sufficiently well represented by a rotating disc or 
whether radial variations of PA and INCL were required. 
All the ROTCUR models are circularly symmetric. 

Subsequently the fits we re used as an input for DiskFit 
llSpekkens fc SellwoodI [ 2 OO 7 I ') in a semi-automated routine, 
in order to avoid a bias the ROTCUR fit outputs were varied 
randomly by up to 5° in PA and inclination and 10 km s“^ 
for the systemic velocity to obtain initial guesses for Disk- 
Fit. Each of the galaxies was then fitted with a rotation-only 
model using DiskFit. When the residual velocity fields for 
these fits showed clear indications of noncircular motions, we 
added bi-symmetric flows or a quadratic warp to the mod¬ 
els. For DiskFit we excluded galaxies that are too small or 
that have too low rotation velocities. The models used in 
this discussion all converged. 

From the initial 25 galaxies we fitted 16 reliably with 
DiskFit and 25 with rotcur. Table [5] shows the basic pa¬ 
rameters of these 25 galaxies as found in the literature. The 
last column indicates wether the galaxy was successfully fit¬ 
ted with DiskFit. When fitting with DiskFit and rotcur 
we used velocity fields derived from the data cubes by fit¬ 
ting a Gauss-Hermite polynomial to each line profile (see 
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Figure 4. Rotation curves used and fitted in the artificial galaxy database. Left panel: dwarf galaxy and intermediate galaxy. Right 
panel: massive galaxy. The black circles and lines indicate the input rotation curves, grey solid lines are the reliable output curves, grey 
dashed lines are from artificial galaxies outside the reliable fitting range (see section [4Tj. 



Database Galaxies 

Variations 

Parameter 

Base 1 

Base 2 

Sample Values 

Beams across the major axis 

16 

16 

4, 6, 7, 8, 9, 10, 11, 12 

Inclination (°) 

60 

50 

10, 20, 30, 40, 85, 88, 90 

PA (°) 

45 

55 

- 

Inner-Outer Dispersion (km s“^) 

8.0-8.0 

8.0-8.0 

13.0-7.5 

Inner-Outer A d {°) 

0-0 

0-0 

0-0, 0-9, 0-3 

Inner-Outer A 0 (°) 

0-0 

0-0 

0-17, 0-3, 0-11 

Inner-Outer scale height (kpc) 

0.3-0.3 

0.2-0.2 

0.2-0.5, 0.5-0.5 

Channel width (km s“^) 

4 

4 

- 

Signal-to-Noise 

8 

8 

2,4,16 

Rotation Curve 

intermediate 

massive 

dwarf 


Table 1. Values for the two base galaxies (columns 2 & 3) and the values by which they are varied (column 4). We change only one 
parameter at a time while all other parameters for the base galaxies are unchanged. 


Figure [5l left panels; only galaxies that are fitted with all 
three codes are displayed). 

4 RESULTS 

4.1 Artificial Galaxies 

The results of the test on artificial galaxies are in Figure 
[6] which shows the difference between input and the best¬ 
fitting values returned by FAT. This difference is calculated 
either as the simple difference for parameters that do not 
change with radius, e.g. the dispersion or systemic velocity, 
or as the error-weighted absolute mean deviation for param¬ 
eters that vary with radius, i.e. PA, INCL, SBR and VROT. 
While the former is expected to oscillate about a value of 
zero across the sample as a whole, the latter is intrinsically 
a positive quantity. The star symbols show artificial galaxies 
that lie outside the range in which FAT is deemed reliable. 


This range is determined in an iterative process by identify¬ 
ing galaxies for which a single fit parameter lies more than 
3 (t from the mean absolute deviation over the whole sample. 
Subsequently these galaxies are excluded when calculating 
the averages and rms. The following section explains where 
FAT fails in more detail for each parameter and what the 
final limits are. 

We express the results in Figure |6] in units of beam 
(30”), channel (4 km s“^), noise and beams across the ma¬ 
jor axis where possible so that they may be scaled to other 
emission line observations. 

4.1.1 Central Coordinates & Position Angle 

We see that for the test galaxies the pipeline recovers the 
input of most parameters to within the estimated errors. 
The central coordinates (Figure panel a and b) are 
accurate to a fifth of a channel in systemic velocity and 
the discrepancy in RA and DEC is in general lower than a 
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Figure 5. Velocity fields of the subset of the LVHIS sample that has been fitted with DiskFit, rotcur and FAT. From left to right each 
row shows the velocity field of the Data as well as the best fitting, FAT, ROTCUR and DiskFit models. The first two velocity fields are 
derived from a Gaussian-Hermite fit. The ROTCUR velocity fields are produced from the fitted parameters with the GIPSY task VELFI and 
the DiskFit velocity fields are output from that code. The black contours start at the minimum velocity in the field and have increments 
of 25 km s“^. The black ellipses show the PA and inclination at the furthest extent of the models. 
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Figure 6. Difference between the artificial galaxy inputs and FAT output. On the y-axis the difference between the model input and 
best fitting FAT values is shown for various parameters. The solid line shows the mean deviation over the whole sample and the dashed 
lines show the la offsets from this average. The errors on the points are calculated as the average of the FAT error for all rings divided 
by \/N, where N is the number of rings in the model. The stars show galaxies outside the reliable range for FAT whereas circles are 
inside this range. The colour coding indicates: blue - dwarf galaxy, black - intermediate galaxy, red - massive galaxy.The panels show 
a) central coordinates difference vs. systemic velocity difference, b) declination difference vs. right ascension difference, c) PA difference 
vs. inclination, d) inclination difference vs. inclination, e) rotational velocity difference vs. inclination, f) velocity along the line of sight 
vs. Inclination, g) inclination difference vs. the number of beams across the minor axis, h) inclination difference vs. maximum diameter, 
i) scale height difference vs. inclination, j) dispersion difference vs. maximum diameter, k) surface brightness difference vs. noise, and 1) 
model extent difference vs. model extent. Panels c), d), e), f), g) and h) show absolute differences (see text). 
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LVHIS Galaxies 


HIPASS Name 

Alternative Name 

a 

J2000 


S 

J2000 


Vsys 
km s~^ 

Wso 

km s~^ 

PA 

O 

Inclination 

o 

Distance 

Mpc 

DiskFit 

converged 

HIPASS 

J0047-25 

NGC 253 

OOh 

47m 

31s 

-25d 

17m 

22s 

243± 2 

407 

52 

83±0 

3.1 

Y 

HIPASS 

J0256-54 

ESO 154-G023 

02h 

56m 

55s 

-54d 

34m 

58s 

574± 2 

122 

39 

35 88 

6.8 

Y 

HIPASS 

J0317-66 

NGG 1313 

03h 

17m 

57s 

-66d 

33m 

30s 

470± 2 

168 


36±7 

4.0 

Y 

HIPASS 

J0320-52* 

NGG 1311 

03h 

20m 

05s 

-52d 

11m 

34s 

568± 5 

80 

40 

79±1 

5.3 

Y 

HIPASS 

J0333-50* 

IC 1959 

03h 

33m 

15s 

-50d 

25m 

17s 

640± 4 

128 

147 

88±4 

8.2 

N 

HIPASS 

J0705-58 

AM 0704-582 

07h 

05m 

18s 

-58d 

31m 

19s 

564± 2 

68 


65I 

4.9 

Y 

HIPASS 

J0731-68 

ESO 59-GOOl 

07h 

31m 

20s 

-68d 

llm 

19s 

530± 3 

82 


20±18 

4.5 

Y 

HIPASS 

J1047-38* 

ESO 318-G013 

lOh 

47m 

39s 

-38d 

51m 

45s 

711± 7 

42 

75 

35 88 


N 

HIPASS 

J1057-48 

ESO 215-G7009 

lOh 

57m 

32s 

-48d 

llm 

02s 

598± 2 

67 

72 

64±27 

5.3 

Y 

HIPASS 

J1219-79* 

IC 3104 

12h 

19m 

04s 

-79d 

42m 

55s 

429± 4 

40 

45 

59±6 

2.6 

N 

HIPASS 

J1305-40* 

GEN 06 

13h 

05m 

02s 

-40d 

06m 

30s 

617± 4 

33 


60^ 

6.1 

N 

HIPASS 

J1305-49 

NGC 4945 

13h 

05m 

24s 

-49d 

29m 

35s 

563± 3 

361 

43 

85±4 

4.1 

Y 

HIPASS 

J1321-36 

NGC 5102 

13h 

21m 

55s 

-36d 

38m 

03s 

468± 2 

200 

48 

70±6 

18.8 

Y 

HIPASS 

J1337-28** 

ESO 444-G084 

13h 

37m 

18s 

-28d 

02m 

17s 

587± 3 

56 


39±4 

5.1 

N 

HIPASS 

J1337-29 

NGC 5236 

13h 

37m 

02s 

-29d 

55m 

03s 

513± 2 

259 


2±27 

6.8 

Y 

HIPASS 

J1337-39* 


13h 

37m 

30s 

-39d 

52m 

56s 

492± 4 

37 


36i 

4.8 

N 

HIPASS 

J1337-42 

NGC 5237 

13h 

37m 

47s 

-42d 

50m 

51s 

361± 4 

77 

128 

35±0 

3.7 

Y 

HIPASS 

J1348-53 

ESO 174-G7001 

13h 

48m 

01s 

-53d 

21m 

31s 

688± 3 

71 

170 

76±11 


Y 

HIPASS 

J1403-41** 

NGC 5408 

14h 

03m 

21s 

-41d 

22m 

26s 

506± 3 

62 

62 

55±8 

4.9 

N 

HIPASS 

J1428-46* 

UKS 1424-460 

14h 

28m 

06s 

-46d 

18m 

32s 

390± 2 

48 


73I 

3.4 

Y 

HIPASS 

J1413-65 

Circinus 

14h 

13m 

27s 

-65d 

18m 

46s 

434± 3 

242 

40 

64±4 

4.2 

Y 

HIPASS 

J1441-62* 


14h 

41m 

37s 

-62d 

44m 

38s 

672± 8 

52 




N 

HIPASS 

J1501-48 

ESO 223-G009 

15h 

01m 

08s 

-48d 

17m 

04s 

588± 2 

61 

135 

44±19 

6.0 

N 

HIPASS 

J1620-60* 

ESO 137-G018 

16h 

20m 

56s 

-60d 

29m 

18s 

605± 3 

139 

30 

73±2 

5.9 

Y 

HIPASS 

J2202-51 

IC 5152 

22h 

02m 

41s 

-51d 

17m 

37s 

122± 2 

84 

100 

51±4 

1.8 

Y 


Table 2. Overview of the 25 LVHIS galaxies selected for testing the code. All positional data is taken from Table 2 in iKoribalski et al.l 
l|20Q4l b Th e distance i s the currently stated average distance in the NASA Extra-galactic D atabase (NED). The P A and axis ratios are 
taken from [haubert j lll982t) . The inclination is calculated fro m the stated axis ratios (see lAaronson et aflllQSOh and the error is the 
difference between the value calculated from the axis ratios in iLaubertsI lll982ll and those provided in NED. If no error is stated this 
means that the values were the same or solely from NED (marked ‘Galaxies with less than eight beams across the major axis in FAT 
fits. “The FAT failed for these galaxies. 


tenth of a beam. 

The PA (Figure [6l panel c) has a mean absolute 
deviation from the input values of 0.23°. It is clear from 
this panel that this deviation is in general much smaller 
than the estimated errors and hence we deem the PA to be 
fitted successfully for all artificial galaxies. 


4-1.2 Rotation and Inclination 

The second row from the top of Figure |6] shows the results for 
the inclination (panel d), rotation (panel e) and the line of 
sight velocity (vrot x sin(i)) (panel f). It is immediately clear 
from the last panel (panel f) that the line-of-sight velocities 
are recovered to within a channel. The trend with inclina¬ 
tion is caused by the weighted average of the deviation: at 
lower inclinations, the range in line-of-sight velocities be¬ 
comes smaller and smaller and hence the deviations become 
smaller. 

The fact that these line-of-sight velocities are well fit¬ 
ted means that the large deviations in rotational velocities 
(panel e) at low inclination (< 20°) are caused by the er¬ 
rors in inclination. This is understandable as an error of 1° 
in inclination immediately leads to an offset of ~ 4 km s“^ 
(~1 channel in this test sample) in rotational velocity for a 
rotation velocity ~ 200 km s“^. The code however amplifies 


these errors, and it is not clear why. 

There are four galaxies with rotation curves that devi¬ 
ate from the actual values by more than 3a when consider¬ 
ing the rotation curve. Clearly the fit is no longer reliable 
at 10° inclination. The outlier at 50° is an artificial galaxy 
with only 4 beams across the major axis. This is also visible 
in Figure 0] where the two dashed rotation curves for the 
massive galaxy represent these galaxies. The fourth outly¬ 
ing galaxy is the massive galaxy at 30° inclination; it has 
a deviation that is a tenth of a channel more than 3 ct from 
the mean deviation (a = 0.5 channel and the mean devia- 
tion=0.7 channel) when the error (0.2 channel) is included. 
It seems that this galaxy falls victim to the success in the 
other galaxies and is only an outlier due to the very small 
scatter in the other fits. Since the galaxies at i =20° and the 
intermediate galaxy at i =30° are fitted well we do not ex¬ 
clude this outlier, hence it appears as a red dot rather than 
a star in Figure |6] From this test we conclude that FAT can 
reliably fit rotation curves in a range from 20°-90° inclina¬ 
tion. 

This lower limit in inclination implies FAT can model 
galaxies that a re more face on than 2D methods can (40°; 
lBegemanlll987l l. It is possible that FAT can fit galaxies at 
even lower inclinations reliably; however, our simulations do 
not probe much below this threshold and therefore we adopt 
20 ° as a conservative lower limit. 
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We also highlight that, while the errors at low inclina¬ 
tions are high compared to what is reachable for higher in¬ 
clinations, they might still be statistically unbiased, as panel 
e) (line-of-sight velocity) suggests. Another test series with 
a significantly larger sample size is required to determine if 
this is indeed the case. 

In inclination (Figure [6l panel d) the typical error is ~ 
1° and the scatter in accepted fits is 0.8°. The two outliers at 
inclinations of 50° and 60° are galaxies with 4 beams across 
the major axis. This is also visible in panel h of Figure [6] 
which shows the difference in inclination plotted against the 
number of beams across the major axis. From this plot it 
is immediately clear that, even in 3D, the inclination esti¬ 
mate is hampered by the number of beams across the major 
axis. There is a clear trend of increasing absolute deviation 
and error increase with a decreasing number of beams across 
the major axis. Even though the deviations are both posi¬ 
tive and negative in the sample as a whole, for these smaller 
galaxies the deviations are all positive. 

This behaviour also affects the rotation curve but not as 
strongly as might be expected. Only in the case of the mas¬ 
sive galaxy with four beams across the major axis does the 
deviation from the input stand out, i.e. the outlier in panel 
e at 50° inclination. However, an inspection of the rotation 
curves in these galaxies shows that in general the shape of 
the inclination and rotation curve are retrieved correctly, al¬ 
beit with an offset. Only at 4 beams across the major axis 
does the fit fail completely, returning an inclination that is 
far offset from the input value. The outliers at low inclina¬ 
tion were already deemed unreliable from their deviations 
in rotational velocity. 

Figure 0] shows that the best fitting rotation curves 
of two massive galaxies have significant deviations at large 
radii. The curve with an outer decrease corresponds to a 
warped artificial galaxy. Here the variation in inclination 
cuts through zero and subsequently increases again beyond 
the cutoff radius. This causes the inclinations in the outer 
radii to be overestimated thus lowering the rotation curve 
amplitude. The curve showing the large outer increase cor¬ 
responds to a massive artificial galaxy with 8 beams across 
the major axis. This behaviour is caused by the fit being ex¬ 
tended beyond the truncation radius, with the consequence 
of an arbitrary fitted value for the rotation velocity. It is 
not clear why this happens in this galaxy but not others, 
but clearly panel 1 in Figure [6] shows that this is the most 
extreme case. 

An individual inspection of the fits showed that at high 
inclination the fits for the massive galaxy are not as good as 
their overall deviations indicate. The galaxies at 88° and 90° 
have inner inclinations which are underestimated by ~3°-5°. 
This does not happen for the intermediate galaxies; there¬ 
fore it is likely related to the shape of the rotation curve. 
This clearly points to difhculties with fitting of inner parts 
of edge-on galaxies which is not unexpected. This behaviour 
should be reflected by an increase in the error towards the 
inner parts of these highly inclined galaxies. However, the 
current version of FAT fails to pick this up. As this does not 
affect any other aspects of the fit and the deviations are still 
acceptable we do not exclude this range. Even so, we stress 
that care should be taken when FAT is applied to highly 
inclined (INCL 88°) galaxies. 

Panel g in Figure 0] shows the difference in inclination 


plotted against the number of beams across the minor axis. 
This panel makes it clear that in these 3D fits the errors are 
truly driven by inclination or the number of beams across 
the major axis and not by the number of beams across the 
minor axis. 


4-1-3 Scale Height 

Panel i of Figure 0] shows that on average the scale height is 
over estimated by roughly a third of the beam. The typical 
distance for this test set of galaxies is 19 Mpc. This 
means that the scale heights of the models are about a 10th 
of the beam. The fact that the problem is worse for the 
massive galaxy, which has a smaller intrinsic scale height, 
points also to the scale height being unresolved. However, 
the trend of increasing deviations towards lower inclination 
points to the scale height being less constrained at lower 
inclinations. This is expected, since the most important 
effect of an increased scale height, a thickening of the 
projected disc along the minor axis, scales with the sine of 
the inclination. 

In order to determine the exact cause of this error a 
larger database more specific to this problem is required. 
A comparison with the un-regularised fits hints, at least at 
low inclination, to the various asymmetries and noise peaks 
being compensated in this parameter as the discrepancies 
between fit and model are much lower in the un-regularised 
cases. Like the inclination, we notice the tendency to stretch 
the galaxy model along the minor axis. We hence conclude 
that scale heights less than 1/lOth of a beam are not fitted 
accurately. 


4-1.4 Dispersion 

The fitted dispersion is in general overestimated by half a 
channel (panel j in Figure [G]). Even though half a channel 
is reasonable as the error one can expect, the systematic 
nature of the offset is remarkable. The trend of increasing 
deviation with decreasing number of beams across the ma¬ 
jor axis points towards a beam smearing effect. This idea is 
strengthened by the fact that the intermediate mass galaxy, 
which has a rotation curve with strong curvature, is more 
affected than the massive galaxy with a constant or approx¬ 
imately linear rotation curve. Since TiRlFiC models a ro¬ 
tation curve by linearly interpolating the rotation velocity 
between ring radii, the average deviation of the true rotation 
curve and the modelled one is much higher for the interme¬ 
diate mass galaxy than for the high mass galaxy. We suspect 
that TiRlFiC compensates for this by increasing the disper¬ 
sion more for the intermediate mass galaxy and less for the 
high mass galaxy. 

FAT does not provide an error on the dispersion but 
we estimated an error on each fit of 0.34 channel from the 
scatter in the fits with 16 beams across the major axis. This 
means that the intermediate galaxy with six beams in diam¬ 
eter has a deviation in dispersion more than 3cr larger than 
the mean absolute deviation in the sample. As this is clearly 
a continuous trend we conclude that FAT becomes unreli¬ 
able for galaxies with less than seven beams across the major 
axis. This is a conservative limit as clearly this depends on 
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the shape of the rotation curve. 

At this stage it is unclear whether the continuous over¬ 
estimation of the dispersion is an artefact of the code, or 
is related to how the artificial galaxies are constructed, or 
whether the small deviations in other parameters are com¬ 
pensated by offsets in dispersion. This latter option does not 
seem to be the case though as the dispersions are similar for 
the regularised and unregularised fits. 

Of additional interest here are the two galaxies where 
the input dispersion is varied from 13.0 km s“^ in the in¬ 
nermost radius to 7.5 km s“^ at the outer most radius. In 
case of the intermediate mass galaxy the code fits a total 
dispersion of 12.9 km s“^ indicating that the fit is driven by 
the central regions of the galaxy. However, for the massive 
galaxy the code fits 10.3 km s“^ which is intermediate to 
the two extremes. This galaxy is also the only galaxy with a 
positive offset because for single values, such as the disper¬ 
sion, only the value in the central ring is considered. This 
result means that from these galaxies it is unclear how the 
code deals with varying dispersions and a larger test set is 
required. 

4-1.5 The Surface Brightness & Signal-to-Noise 

The surface brightness is captured well in general as shown 
in panel k of Figure [S] The average deviation is ~ 4%. This 
panel shows the difference between the flux in the input 
model and FAT output, as a fraction of the total flux in the 
input model, against the number of beams across the major 
axis. The total flux was measured in the original model 
cubes without added noise, but with the correct resolution, 
in voxels where the flux exceeded the subsequently added 
noise. The error is the formal error expected from the noise 
and number of added resolution elements. 

It is clear from Figure [6] panel k that the number 
of beams across the major axis severely affects the fitted 
flux. In this case galaxies that lie more than 3a from the 
mean are those with four beams across and the massive 
galaxy with six beams across. These deviations appear 
to be connected to those in dispersion. However, a larger 
deviation in dispersion seems to lead to a lower deviation 
in total flux, or vice versa. At this point it is unclear how 
these two parameters correlate with each other although 
the clear increase with number of beams across the major 
axis points towards a resolution effect. 

FAT underestimates the total flux in all models. This 
is partially due to the models being constructed at lower 
resolution than the artificial galaxies. This means that the 
FAT models are unable to capture the high flux peaks in the 
artificial galaxies. From this dataset it is unclear whether 
this fully explains the offset. However, there is also a small 
trend with inclination (not shown) where the fits at lower 
inclinations capture the total flux much better than those 
at high inclination. This confirms the idea that the offset 
is mostly caused by the peak flux being underestimated as 
this ring of high flux would cover a larger sky area in lower 
inclination galaxies and would thus be easier to fit. 

Besides the trends with the number of beams across 
the major axis and inclination the total flux also shows a 
weak trend with signal-to-noise. The two points with the 
largest deviation are the galaxies with a signal-to-noise of 
two. This is expected as the lower signal-to-noise would 


make it harder to determine the correct flux. The lower 
signal-to-noise hides more of the flux in the noise and hence 
a model with lower total flux is fitted to the input. However, 
the deviation of these galaxies is still acceptable. 


4.1.6 Extent of the models 

The last panel in Figure (panel 1) shows the difference in 
adopted outer radius vs. input outer radius. FAT estimates 
the size of the galaxies in this test sample quite accurately. 
Since the edges in the artificial galaxies are created by 
decreasing the scale length by a factor of 20 at the cutoff 
radius, FAT can be used to detect these edges in real 
galaxies if they exist. However, as at the given edge the 
emission does not immediately disappear but continues into 
the noise, FAT sometimes overestimates the outer radius 
by 1 beam. This is not unexpected. There appears to be 
a trend in which the overestimation increases for galaxies 
with fewer beams across the major axis. 

From these tests we conclude that FAT fits galaxies well 
in the inclination range 20°-90° as long as the errors pro¬ 
duced by the pipeline are taken into account. At an inclina¬ 
tion of 10° the errors in rotation velocity become so large 
that they render the derived rotation curve meaningless. 
However, considering the derived errors, the true rotation 
curve is still reproduced accurately in a purely statistical 
sense. The same is true for the inclination. The minimum 
size that can be reliably fitted without any priors is a fi¬ 
nal model with 7 beams across the major axis. In galaxies 
that are less resolved the dispersion or surface brightness is 
seriously misfitted, most likely to counterbalance the fail¬ 
ure of the current TiRlFiC implementation to implement 
curvature, between rings, in a small scale model. However, 
because FAT overestimates the size of the galaxy to a larger 
extent when it is smaller, the models fitted to artificial galax¬ 
ies with six beams across can have 8 beams across the major 
axis in the output. As we want to determine reliability from 
the fitted models we conclude that without additional input 
only FAT models with eight beams or more across the ma¬ 
jor axis and between 20°-90° inclination are reliable. This is 
a very conservative limit because the tests here show that 
fits to smaller galaxies (up to six or seven beams across the 
major axis) can be accurate. 

It is clear from Figure [6] that FAT is slightly better at 
fitting the massive galaxy (red) than the intermediate mass 
galaxy (black). This affects the fitted dispersion the most 
but is also visible in other parameters. Further tests are re¬ 
quired to see whether this is due to the shape of the rotation 
curve or because of the higher rotational velocities. The two 
dwarf galaxies simulated here unfortunately do not shed any 
further light on this matter. 

4.2 LVHIS Sample 

TiRiFiC and FAT flag 25 of the 25 LVHIS galaxies as 
successfully fitted. A visual comparison between the FAT 
models and the data cubes confirms that the intensity dis¬ 
tributions of all models, barring two (HIPASS J1337-28 and 
HIPASS J1403-41), correspond to the data. For HIPASS 
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Figure 7. Differences between the final fits with FAT, ROTCUR and DiskFit. Left hand side of each panel: the difference between ROTCUR 
and DiskFit, i.e. rotcur - DiskFit, Right hand side of each panel: blue: the difference between ROTCUR and FAT, red: the difference 
between DiskFit and FAT. The solid grey lines show the weighted average and the dashed lines show the Icr deviation from this average 
in each plot. Stars correspond to final fits outside the previously determined reliable range (see Section f4.1ll and circles to ones inside these 
limits. The errors on the points are calculated as the average of the error for all rings divided by \/N with N the number of rings in the 
model. The panels show a) Central coordinate difference vs. Systemic velocity difference, b) Right Ascension difference vs. Declination 
difference, c) PA difference vs. Inclination, d) Inclination difference vs. Inclination, e) Rotational velocity difference vs. Inclination, f) 
Velocity along the line of sight vs. Inclination, g) Rotational velocity difference vs. Model extent, h) Model extent difference vs. Model 
extent. Panels c), d), e), f) and g) show absolute differences (see text). 
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Parameter 

MeanRC-DF 

f^RC-DF 

MeanRC-FAT 

CTRC-FAT 

MeanoF—FAT 

ctdf-fat 

RA (beam) 

-0.24 

0.19 

0.09 

0.79 

0.22 

0.82 

DEC (arcsec) 

-0.24 

0.48 

-0.09 

0.27 

0.13 

0.40 

Vsys (km s“l) 

-0.02 

1.34 

0.13 

0.92 

0.17 

0.83 

PA (°) 

2.75 

5.79 

7.41 

4.53 

5.81 

4.18 

Inclination (°) 

1.72 

3.00 

10.47 

5.91 

11.71 

6.59 

Vrot (km s“l) 

0.99 

0.65 

3.17 

2.41 

3.66 

2.43 

VrotXsin(i) (km s“^) 

0.46 

0.64 

0.56 

0.39 

0.70 

0.58 

Max. Diameter (beam) 

0.96 

2.85 

-3.11 

3.56 

-5.23 

6.06 


Table 3. The mean absolute deviation and rms scatter of the difference between the classical ROTCUR and DiskFit fits, the ROTCUR and 
FAT fits and the DiskFit and FAT fits for galaxies in the reliable range. For the first and last 2 columns this is for 11 galaxies and the 
middle 2 columns is a sample of 22 galaxies. 


J1403-41, the kinematics are co mplex and not well repro - 
duced by a tilted-ring model (see Ivan Evmeren et al. IH), 
while for HIPASS J1337-28 imaging artefacts misguide the 
fitting procedure. These two galaxies are marked with ** in 
Tableland excluded from further analysis in this section. 
Optimally, FAT should be able to provide the metrics to 
exclude such cases. This task will be the subject of future 
research. 

In the following we compare the parameters obtained 
from FAT to those obtained from ROTCUR and DiskFit 
fits as well as to the literature in order to get an idea of the 
errors and differences between the codes. 

An overview of the results of the fits by FAT for the 
sample selected from LVHIS is given in Figure [3 and a 
subsample of velocity fields for the data and the three 
codes is shown in Figure [S] The panels in Figure [7] show 
the difference between the ROTCUR and DiskFit fits (black 
symbols, left side of a panel), the difference between the 
ROTCUR and FAT fits (blue symbols, right side of a panel), 
and the difference between the DiskFit and FAT fits (red 
symbols, right side of a panel) on the y-axis plotted against 
various parameters. As for the artificial galaxies, in cases 
where the values may vary with radius we calculate an 
absolute weighted mean deviation of the ring differences and 
simply take the difference between values that do not vary 
with radius. As the fits were sometimes performed using 
different ring sizes, we performed a linear interpolation of 
the ROTCUR or DiskFit fit to the FAT radii when needed. 

The errors on the points were calculated in the follow¬ 
ing manner. For parameters that have errors the errors for 
for both codes were added in quadrature and subsequently 
averaged and divided by the square root of the number 
of rings. For parameters that were provided without error 
no error is shown in Figure [T] Table [3] shows the average 
deviation of the total sample and the rms of the sample that 
falls within the previously determined reliable FAT fitting 
range. These values are also indicated as the solid grey 
lines (average) and the dashed lines (Icr deviation from the 
average) in Figure [3 Galaxies that fall outside the range 
determined in section Hr] are marked by star symbols. 

It is important to note that the FAT, ROTCUR and 
DiskFit models do not always have the same number of 
degrees of freedom. For example, sometimes the researcher 
chooses not to include a warp in the ROTCUR model where 
FAT decided to fit a warp; also, it can happen that in 
the DiskFit models streaming motions are included. We 
compare these final models because we want to investigate 


how these different choices and methods affect the fits. 
Additionally, we want to compare these methods at their 
best and minimise the error on the fits. In the following 
sections we will discuss the fitted parameters in more detail. 

Applying the conservative limits derived in section 
[4T] thirteen galaxies fitted with DiskFit and fourteen 
galaxies fitted with ROTCUR can be fitted reliably with 
FAT. We include the unreliable fits in this analysis for 
several reasons. They serve as an illustration of why such 
galaxies should not be fitted and show that the regularity 
of a velocity field is no guarantee for a reliable fit. It is 
also interesting to see where the comparison between 3D 
and 2D fits starts to de viate as a con firm ation of the relia - 
bility limits proposed bv iBosm ili^) an d lBeeem^ lll987li . 


4-2.1 Central Coordinates 

First we will look at the central coordinates (Figure [3 
panels a and b). We find that FAT finds centres comparable 
to ROTCUR and DiskFit. The grey circle in Figure [3 (panel 
b) for the A RA and A DEC plot delineates the ellipse 
with (TARA as the horizontal axis and (tadbc as the vertical 
axis (see Table [3|. These are on average half of the FWHM 
for the observations in this sample. Figure [3 illustrates that 
the centres fitted by the 3D branch agree better with the 
DiskFit fits than with the rotcur fits, even better than 
DiskFit agrees with rotcur. 

The single outlier in RA is HIPASS J0317-66. It turns 
out that in this case it is difficult determine a single centre. 
The centre found in HIPASS (see table [2]) disagrees with 
the DiskFit, rotcur and FAT fit (d = -66°29'37.8", 
S = -66°30'4.0", S = -66°29'52.1", a = 3‘'18“17.6D, 
a = 3** 18™ 17.13", a = 3‘'18™9.47", respectively). The 
optical centre agrees with the centre found in the ROTCUR 
fit but a visual inspection of the data cube favours the 
centre found by FAT. As this is a large and complex galaxy 
(panel b in Figure [5]), and a more detailed inspection of the 
data does not favour one centre over the other, we ignore 
this mismatch in our interpretation. 

For the systemic velocity (panel a) we see that although 
the FAT fits for some galaxies differ by more than a channel 
width from the 2D results, the scatter is comparable to 
that between ROTCUR and DiskFit. We therefore conclude 
that FAT successfully fits the centres of these galaxies. 
From this we estimate that a typical error on the systemic 
velocity is about a channel in real observations, and that 
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that on the spatial coordinates is about a fifth of the beam. 


4-2.2 Position Angle & Inclination 

Similar to the central coordinates we see that FAT provides 
a good fit for the PA compared to ROTCUR and DiskFit 
(Figure [71 panel c). It is notable however that the scatter 
for this sample is ~ 4° and the mean absolute deviation 
~ 6.6°. This is due to a combination of errors and different 
fitting methods, as the position angle is allowed to change 
in the centre for the ROTCUR fits but not for the FAT or 
the DiskFit fits. The largest outlier in this plot is HIPASS 
J1413-65. For this galaxy the ROTCUR fit to the PA does 
not vary whereas the FAT model shows a variation in PA 
with several directional reversals, i.e. the PA increases, 
decreases and subsequently increases again. This leads to 
the large absolute mean deviation. A warp does seem to be 
favoured by the data (see Figure O panel g) and the 3D 
model traces the emission in the data cube well. 

It is clear from Figure [7] (panel d) that the biggest 
discrepancies can be found when comparing the inclina¬ 
tions. However, the deviations in the line-of-sight velocities, 
i.e Vrot X sin(i) (panel f in Figure [7]), are similar in most 
models and the scatter in their difference is half a channel 
(cTobs^ 0.5 channel). This makes it clear that FAT, DiskFit 
and ROTCUR all solve for the degeneracy between rotational 
velocity and inclination in a different way. The outlier is 
a large edge-on galaxy (HIPASS J1305-49, Figure [5] panel 
p) for which the 2D codes fit an inclination 10° lower than 
does FAT. However, as this happens at high inclination the 
rotation curve is hardly affected leading to a large offset in 
the line-of-sight velocity. 

From panel d in Figure [7] it can be seen that DiskFit 
and ROTCUR overall agree on inclination whereas FAT 
clearly deviates at all inclinations. Given the results in 
section 14.11 we cannot exclude systematics errors in FAT 
for galaxies smaller than 8 beams in diameter even in 
idealised tests. Hence, we will focus on the galaxies in the 
range reliably fitted by FAT in what follows. 

Relative to 2D methods FAT stretches the inclination 
range that is fitted. At low inclination FAT finds inclina¬ 
tions lower than the 2D methods and at high inclination 
ones that are higher (see Figure (SI). Note that only the 
velocity fields for the galaxies fitted with DiskFit are 
shown in Figure [5] as the other galaxies are more complex 
and the interpretation of the velocity field is not straight 
forward. In these velocity fields FAT outer rings do not 
always correspond to the actual map because the fitting 
is done on an unmasked cube and can go below the noise 
in certain cases. At the low-end of the inclination range 
considered here, the ellipses and velocity fields show a lower 
inclination for the FAT fits as they are more circular. It is 
hard to differentiate by eye between the different models 
and therefore we also present the residual velocity fields in 
Figure (8] 

These residual fields emphasise that the codes solve for 
the degeneracy between inclination and rotation in different 
ways. Even though the inclinations differ significantly, the 
residuals are similar. However, is the 3D method doing 
better? In order to address that question we will examine 
an individual galaxy in more detail. 


Figure [9] shows the velocity field, an integrated mo¬ 
ment map and an infrared image of HIPASS J0256-54. The 
difference in inclination between the DiskFit and rotcur 
fits with FAT is 17° and 16°, respectively. A closer look 
at the inclinations (Figure HOI) of the models shows that 
the DiskFit and rotcur fits are fiat at ~ 67° inclination 
whereas the 3D model is warped and varies from 89° in 
the centre to ~70° in the outer parts. It is clear from the 
optical image and the integrated moment map in Figure [5] 
that this galaxy is close to edge-on in its inner parts, indeed 
from the optical axis ratios an inclination ^ 88° (see Table 
mi is expected. From the residual velocity fields it seems 
that the FAT fit has somewhat lower residuals (see Figure 
(8] panel o) but not significantly so. 

Figure [TT] shows the difference between the literature 
values for the inclination and PA, as stated in Table [2] and 
all of the models. To ensure an appropriate comparison, the 
PA and inclination of the Hi observations were determined 
at the optical major axis radius found in the literature. 

From the left panel in Figure[Tl]it becomes immediately 
obvious that to determine a reliable PA at low inclination, 
kinematical information is crucial as at the low inclinations 
there is a large difference between the optical PA and the 
PAs determined from the fits, regardless of the method. 
Below an inclination of 40° this diffe rence is due to the 
inaccuracy of the optical PA estimates llLaubertslll982ll . 

When we turn to the right hand panel of Figure [TT] it 
appears that at the high inclinations FAT performs slightly 
better than the 2D methods. However, this is based on 
only three galaxies (HIPASS J0047-25, HIPASS J0256-54, 
HIPASS J1305-4) and the difference is not significant when 
compared to the overall scatter. An individual inspection of 
HIPASS J0047-25 and HIPASS J1305-4 (panels m and p in 
Figure [5]) shows that the differences are created in a similar 
fashion as previously described for HIPASS J0256-54. In the 
case of HIPASS J0047-25 the optical major axis is actually 
larger than the models and hence the warping in FAT leads 
to a larger discrepancy compared to the 2D methods at the 
radius of the optical major axis. Overall however the fits 
are similar. 

Figure [TT] shows that in the intermediate inclinations 
(40° < i < 70°) the 3D fits, on average, agree with the 
optical values. This is also true for the 2D fits, however from 
Figure [7[ it is clear that there is an average difference of 
~15° between the FAT and 2D inclinations. Once more this 
is caused by FAT fitting warped distributions where the 
2D methods fit fiat disks. One example of this is HIPASS 
J0705-58 (Figure [5] panel h). From this velocity field it is 
clear that beam smearing, warps and a larger extent of the 
model can easily explain these differences. It appears the 
methods do equally well in this range when compared to 
the optical information available. 

At the low end of the low inclination range (< 40°) 
it is clear that the 3D method prefers lower inclinations 
compared to the 2D methods (Figure[7| panel d) and that in 
general both the 2D and 3D methods find inclinations lower 
than the literature values (Figure [TT]). Of the reliably fitted 
galaxies, four agree well between the different methods 
but one does not. This galaxy, HIPASS J1348-53, has 
a difference of more than 20° between the 2D and 3D 
methods. It clearly has an optically determined inclination 
that is much higher than suggested by the velocity field 
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Figure 8. Velocity field residuals of the subset of the LVHIS sample that has been fitted with DiskFit, rotcur and FAT. From left to 
right each row shows the velocity field of the Data, the residuals of FAT, the residuals of ROTCUR, the residuals of DiskFit. The black 
contours always correspond to the contours of the data velocity field. They start at the minimum velocity in the velocity field and are 
drawn at intervals of 25 km s~^. The grey solid line outlines the extent of the data velocity field. 


© 2013 RAS, MNRAS OOO.flM^ 



































































































































Automated Kinematic Modelling in Large Hi Surveys 

-4.0E+18 Njn(cni'^ 3.2E+21 506.00 (kms‘‘) 634.00 


17 


-54-40' 


8 

Q 


-54-40 - 



30.0’ 57’°00.0‘ 30.0' 02''56"00.0‘ 30.0' 57'”00.0' 30.0’ 02'‘56'°00.0’ 30.0’ 57’°00.0‘ 30.0' 02'‘56‘°00.0‘ 


RAa2000^ 
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in the integrated moment map are at 6, 12, 24 .... etc. xlO^® cm~^ and for the velocity fields at 506, 531, 556, 581, 606 km s“^. This 
galaxy is in also shown in panel o in Figures [5l and [8l 


Fitting Code 

PA Mean Abs. Dev. 

PA cr 

Inc. Mean Abs. Dev 

Inc. cr 

No. Gal. (Pa, Inc.) 

ROTCUR 

4 

40 

8 

24 

10,14 

DiskFit 

2 

30 

7 

22 

9,13 

FAT 

7 

38 

4 

24 

10,14 


Table 4. The mean absolute deviation and rms scatter of the difference between literature values and value found from the kinematic 
fitting for the different codes. The values of the model are extracted at the radius corresponding to the major axis radius found in the 
literature. 


(76°, FigureO panel c) making it impossible to differentiate 
between the methods. 


4-2.3 Rotation curves 

As tilted-ring modelling is pre-dominantly used for the 
extraction of rotation curves, we want to focus on how well 
the rotation curve is retrieved for each galaxy. Figure [12] 
shows the rotation curves for the 23 galaxies in our sample 
fitted by FAT. Here the black lines show the rotation 
curves derived with FAT, the blue dashed lines derived 
with ROTCUR and the red dot-dashed lines derived with 
DiskFit. We see that in general, the FAT rotation curve is 
consistent with those determined by the 2D methods albeit 


sometimes with an offset. Comparing these plots to panel d 
in Figure[7]it is evident that these offsets are a result of the 
different inclinations found by the 2D methods. 

In general the fitted rotation curves seem reasonable 
within their errors. A possible difference between the 2D 
methods and FAT is that in several cases the rotation curve 
found with the 3D method rises faster in the inner point. 
However, this also occurs in the inner most point of the 
rotation curves of the intermediate mass artificial galaxies 
(see Figure|4|| and is therefore, at least partially, an artefact 
of the code. Again, TiRlFiC interpolates linearly between 
the single ring radii. Any quantity with a high curvature is 
hence ill-fitted. This results in an overshooting at the fitted 
nodes. 

It is remarkable that also in this sample large devia- 
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Figure 10. HIPASS J0256-54: Parameters for the Rotation curve, 
PA and Inclination. Black + Blue: FAT both sides independently, 
red: ROTCUR green: DiskFit. 

tions start to occur when there are less than eight beams 
across the major axis. This would imply that the methods 
for retrieving rotation curves from veloc ity fiel d s hav e 
significantly improved since the work of iBosmal lll978ll . 
However, a more detailed analysis than provided in this 
paper is required to confirm this. 


5 DISCUSSION 

We have shown that the code presented in this paper can, 
for the first tim43i successfully fit tilted-ring models to Hi 
data cubes in a fully automated manner over a large range of 
projected and physical sizes. It works well for dwarf galax¬ 
ies, large spirals and galaxies in between as long as they are 
resolved by at least eight beams across the major axis and in 
the inclination range 20°-90°. In order to achieve a success¬ 
ful fit all the code requires is the data cube of the galaxy. 
It will then g o on to determine initial estimates (through a 
run of SoFiA. ISerra et ^l2015h . the extent of the model, fit, 
adjust and smooth tilted-ring models without any input or 
adjustments from the user. 

The limits stated here pertain to fits without any pri¬ 
ors on the kinematical parameters - in particular, it is likely 
that galaxies with even fewer resolution elements across the 
major axis could be reliably fit if the inclination were fixed 

® While this paper was under review Bbarolo, another code c a- 
pable of such fits, was released dPi Teodoro fc Fraternalill201fl) . 


to the optical value (e.g. iDi Teodoro &: FraternalilboiSl i. 

The aforementioned limits on the size of a galaxy’s ma¬ 
jor axis relative to the beam and on its inclination in order 
to obtain a reliable FAT model are not hard, but will depend 
on other properties such as, for example, S/N and the shape 
of the rotation curve. Additionally it must be emphasised 
that these are the limits for FAT which is a fully automated 
code. This does not mean that 3D methods cannot fit out¬ 
side these limits. However, in the case of low inclination it is 
very hard to see how one would separate incorrect fits from 
correct fits as the minimum in the landscape becomes 
increasingly more shallow. This is easily seen in Figure 1131 
where the left panel shows the position velocity (PV) di¬ 
agram for the massive artificial galaxy at 10 degrees incli¬ 
nation (black contours, greyscale), the same galaxy without 
noise (blue contours) and the fit (red contours) from FAT. 
Even though the fitted rotation curve is, on average, almost 
20 channels off, the modest deviation in inclination results 
in similar observed velocities. 

On the other hand, for the amount of beams across 
the major axis there are differences. The right hand panel of 
Figure [13] shows a PV-diagram along the major axis of the 
massive artificial galaxy with four beams across the major 
axis. The contours are the same as in the left panel. Even 
though the error in the rotation curve is much smaller for 
this galaxy, clearly the fit deviates from the input artifi¬ 
cial galaxy without noise. However, this difference is hidden 
by the noise in this artificial galaxy. This implies that this 
lower limit for number of resolution elements across the ma¬ 
jor axis can be reduced with improved fitting techniques or 
by increasing the signal to noise. Additionally the require¬ 
ment of eight beams across the major axis is overly strict 
as the models fitted to a galaxy with seven beams across 
the major axis end up with the same extent. This means 
that by merely improving the determination of the extent of 
the model the code can be improved to fitting galaxies with 
merely seven beams across the major axis. Furthermore, the 
choice of different interpolation methods to determine the 
inter-ring behaviour of parameters is also likely to further 
reduce the number of beams across the major axis required 
for a reliable fit. 

The requirement of eight beams across the major axis 
can probably also be pushed down by careful manual adjust¬ 
ment of the model, the introduction of priors (e.g. an optical 
inclination) or a parametrisation of the galaxy model. As in 
the last case the degrees of freedom would be reduced this 
would make a more reliable fit but also the information that 
can be retrieved would be severely reduced, e.g. if one would 
use an analytical function for the shape of the rotation curve 
the maximum velocity could be retrieved and possibly the 
steepness of the rising part but investigations into rotation 
curve shapes would be ruled out. However, such methods 
and limits require further investigation. 

For the specific future of FAT one could easily imag¬ 
ine replacing the first fit of a fiat disc with a parameterised 
model in order to gain speed and accuracy. For galaxies with 
less than eight beams across the major axis this would also 
be the final fit and for better resolved galaxies FAT would 
go on to fit an unconstrained model as it does in the current 
version. A concern here would be how much the initial pa¬ 
rameterised model would steer the final results to this first 
model. 
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Figure 11. Comparison between the PA and inclination of the optical values and the different fits, black - FAT, blue - ROTCUR, red - 
DiskFit. Stars are galaxies where the final FAT model is in a part of parameter space that is considered unreliable. The grey dash line 
denote the scatter on the deviations between literature and the ROTCUR interactive fit. The optical values come from Table [2] 


The FAT can deal with extreme cases of warps. Figure 
M shows the velocity fields of the most extremely warped 
galaxy in the LVHIS sample, HIPASS J1501-48. This figure 
demonstrates the strength of tilted ring fitting in analysing 
Hi observations with complex but regular structures. This 
galaxy also illustrates the problem of retrieving rotation 
curves in low inclination galaxies. Even though both the 2D 
and 3D methods provide good fits to the data the derived 
rotation curves (shown in Figure fT^ panel s) are discrepant. 
This is due to the fact that in the ROTCUR model the incli¬ 
nation is kept flat whereas FAT finds a modest warp. 

FAT fails to accurately reproduce an observed data 
cube in two cases, despite the fact that the code flags the fit 
as being successful. Currently the code does not check the 
models for physical consistency and this could be a way to 
identify cases that need to be inspected. 

When comparing the retrieved inclinations of the fits to 
optical inclinations it is unclear whether FAT retrieves the 
inclination better than the 2D methods (Figure fTTI) at high 
inclinations. In two cases it does, but in one case it does not. 
However, this latter case (HIPASS J0047-25) is a special one 
where the optical disc extends beyond the Hi layer. Other 
galaxies in our sample are either too small or complex and 
therefore this result is only tentative. 

In order to get a better estimate of when to use the 3D 
method and when to use the 2D method in the WALLABY 
kinematics pipeline it will be necessary to determine which 
of the two methods is more suitable for fitting observations 
of barely resolved galaxies, i.e. galaxies with less than eight 
beams across the major axis. The easiest way to investigate 
this is to construct a database of artificial galaxies, larger 
than the one presented here, that can be fitted with the two 
branches of the pipeline. In this manner we know the input 
that we are trying to retrieve and we can accurately deter¬ 
mine which method performs better. Such an investigation 
will also provide a much better idea of the errors on the 
final fits and whether there are other areas of the parame¬ 


ter space where one method performs much better than the 
other. Unfortunately, a simple grid with some variations in 
parameter space leads quickly to a set of thousands of galax¬ 
ies. The current version of FAT is not quick enough to deal 
with this in any reasonable time frame. A new version which 
can be run on high performance clusters is in development 
and will be tested on such a dataset. 

From the tests on artificial galaxies we can get a lower 
limit of the errors in the FAT fits. If we take the rms of the 
difference between the input and output we get a minimal 
error (column 1, Table [5]). As the artificial database lacks 
significant asymmetries and a large number of galaxies this 
scatter can be seen as the minimal error that the code will 
always provide. These are a gross underestimate of the real 
errors. A crude upper limit can be estimated from the scat¬ 
ter in the differences between the 2D methods and FAT 
(column 2, Table [5]). This would be a conservative estimate 
of the error as it includes galaxies in a range where 2D ap¬ 
proaches are unsuitable. In order to check this reasoning we 
compare these values to errors provided on tilted ring fits in 
the literature. 

Recently, several individual galaxies have been mod- 
elled with TiRlFiC ( [C^ntile et al’l l2013l : iKamohuis et al.l 
l2013l : Ide Blok et al.l l2014l L In these models the errors are 
estimated by varying each parameter individually until a 
comparison, by eye, between the data and the model clearly 
shows discrepancies. The errors from these publications are 
shown in the last three columns of Table [5] These confirm 
that the scatter in the artificial galaxies is a lower limit on 
the error, whereas the scatter in the LVHIS sample agrees 
with these conservative errors. The exception to this is the 
large scatter in inclination and PA that is found in the 
LVHIS sample. This is due to the fact that the major dif¬ 
ferences between the codes stem from how they deal with 
PA and inclination. DiskFit always fits flat inner discs, 
whereas FAT and ROTCUR use various degrees of warping. 
Additionally, DiskFit and rotcur fit the whole galaxy at 
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Figure 12. Rotation curves for all the galaxies in the LVHIS Sample. The plots are ordered according to their FAT inclination, with 
the lowest inclination in the bottom left. Black lines: FAT, Red Lines:DlSKFlT, Blue lines:ROTCUR. The error bars are the formal errors 
derived by DiskFit and ROTCUR for the red and blue lines respectively. For FAT the error is an empirical estimate derived from the 
regularisation process (see § 12.21 Stars indicate curves outside the reliably fitted range. 


once whereas FAT divides the galaxies in two along the mi¬ 
nor axis. This does not happen for the rotation curve or the 
central coordinates because similar symmetry requirements 
are applied to these parameters. Therefore, the scatter on 
the PA and inclination is likely to be an overestimation of the 
typical error on these parameters. We conclude, that when 
similar fitting methods are applied, the scatter between the 
2D methods and FAT gives a good estimate of the typical 
error on a FAT fit. Future investigations will aim at an in¬ 
dependent determination of the errors, via bootstrapping or 
Markov Chain Monte Carlo (MCMC) methods. 


6 SUMMARY AND CONCLUSIONS 

In this paper we present a fully automated code for fitting 
tilted-ring models (the Fully Automated Tilted Ring Fitting 
Code or FAT) to Hi d ata cubes. This co de is an IDL wrap- 
per around TiRlFiC lljozsa et al ][2ooa) that uses SoFiA 
jSerra et al.ll2013 1 to obtain initial parameter estimates. It 
is intended for use as the 3D branch of a secondary data 
product pipeline for the ASKAP all sky survey WALLABY 
(IKoribalskil 1201 and its northern counter part WNSHS. 
Besides the 3D code itself we also present the concept for 
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Figure 13. The PV Diagrams for the artificial massive galaxy at 10 degrees inclination (left) and the massive galaxy with 4 beams 
across the major axis (right). Contours are at 2.5, 5, 10, 20 x 10^^ cm~^. Blue contours are the input model, Black the input model 
plus noise and the red contours the final fit. 
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Figure 14. The velocity fields for the galaxy HIPASS J1501-48. This galaxy is the most extremely warped galaxy in the sample. From 
left to right it shows the data, FAT, ROTCUR. Black contours appear at 10 km s~^ intervals. 


Parameter 

Artificial Database 

LVHIS Sample 

de Blok et al. ('2014) 

Kamohuis et al. ('2013) 

Gentile et al. f2013) 

(using rotcur) 



No. Galaxies 

46 

14 

1 

1 

1 

RA (beam) 

0.02 

0.8 


0.22 


DEC (beam) 

0.01 

0.35 


0.31 


■Vsys (channel) 

0.09 

0.9 

1.1 

1 


PA (°) 

0.3 

4.4 

0.6 

1 

2 

Inclination (°) 

0.8 

6.3 

1.5-2.5 

1.3 

2 

Vrot (channel) 

0.5 

2.4 

2.5 

1.2 

1.1 


Table 5. Errors on the tilted ring fit as estimated from the artificial database, the 2D to 3D comparison and in recent literature. 
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this pipeline of which t he 2D branch w ill be presented in 
an accompanying paper doir et al ]|2^. 

In order to test the code we have created a set of 
52 artificial galaxies and sel ected 25 galaxie s from the 
nearby galaxy survey LVHIS dKoribalskil l201(t) . Test runs 
on artificial galaxies indicate an inclination between 20° and 
90° and an extent of eight beams across the major axis as 
a highly conservative lower limit to produce an acceptable 
fit. Further investigations will likely push down these 
limits by a large factor. Even now, galaxies at the lowest 
inclinations are rejected only because the large errors in the 
rotation velocity render the fits unusable for rotation curve 
analyses (but not, e.g. for studies of the surface brightness 
distribution). 

Out of a sample of 25 LVHIS galaxies 23 are fitted suc¬ 
cessfully in the sense that their models trace the observed 
emission when comparing data cubes to model cubes. The 
failed fits are due to residual artefacts in the observations 
and complex emission patterns of the Hi. However, the 
current version of the code fails to flag the unsuccessful fits 
and they need to be identified manually. As the models 
for these galaxies differ significantly from the data this 
could still be achieved without the need for any additional 
observations. 

From the fitting of 23 galaxies in the LVHIS sample 
we show that 2D fitting techniques underestimate the in¬ 
clinations of galaxies at high inclination. As in our current 
sample there are only three galaxies at high inclination 
that are fitted with the 2D methods as well as FAT it is 
unclear whether this is a systematic effect. Additionally, as 
this happens at high inclinations where the change of sin(i) 
is small as a function of i this does not significantly affect 
the amplitude of the derived rotation curve. 

FAT retrieves several parameters related to the bright¬ 
ness distribution such as the scale height and brightness 
profile. When the scale height is significantly smaller than 
the beam this results in an overestimation of the scale 
height, and we conclude that hence the scale height cannot 
be estimated using FAT for our test sample. The current 
tests are not sufficient to identify what happens when the 
scale height is resolved. 

Errors on the individual parameters remain ill under¬ 
stood. However, a comparison with recent modelling work 
shows that the scatter in the sample of LVHIS galaxies is a 
good indication of realistic errors. In order to investigate the 
open questions on the code’s performance, a run on a large 
artificial database of galaxies, where the exact structure 
of the galaxies is known, is required. Unfortunately, the 
current code is not fast enough to deal with the requisite 
number of fits. Hence, the next step is to adapt the code to 
run on high performance computing facilities. 
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